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Abstract 

We provide a unified picture of the domain wall spectrum in supersymmetric 
QCD with N c colors and Nf flavors of quarks in the (anti-) fundamental representa- 
tion. Within the framework of the Veneziano-Yankielowicz-Taylor effective Lagran- 
gian, we consider domain walls connecting chiral symmetry breaking vacua, and we 
take the quark masses to be degenerate. For Nf/N c < 1/2, there is one BPS sat- 
urated domain wall for any value of the quark mass m. For 1/2 < Nj/N c < 1 
there are two critical masses, and m**, which depend on the number of colors 
and flavors only through the ratio Nf/N c ; if m < m*, there are two BPS walls, if 
< m < m„ there is one non-BPS wall, and if m > m„ there is no domain 
wall. We numerically determine and Tri** as a function of Nf/N c , and we find 
that m„ approaches a constant value in the limit that this ratio goes to one. 
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Introduction 



In recent years, several groups have studied domain walls in supersymmetric quantum 
chromodynamics (SQCD) in the framework of the Veneziano-Yankielowicz-Taylor (VYT) 
effective Lagrangian [jl]]. This Lagrangian is not an effective Lagrangian in the usual 
Wilsonian sense, but it it is expected to describe the vacuum structure of SQCD. If the 
masses of the quark matter fields are degenerate, the VYT effective Lagrangian takes the 
form of a Wess-Zumino model in terms of the chiral superfields $, representing the glue- 
ball superfield, and X, representing the quark matter. The effective Lagrangian has N c 
chiral symmetry breaking vacua and an additional chiral symmetry conserving vacuum. 
The spectrum of the domain walls connecting two adjacent chiral symmetry breaking 
vacua was found to depend on the number of colors N c , the number of flavors Nf, and the 
common mass of the quarks m. In Ref.[|2|], domain walls were studied in several models 
with Nf = N c — 1. In the cases of N c = 2, 3 and 4, it was found that for small quark mass 
m, there are two Bogomolnyi-Prasad-Sommerfield (BPS) saturated domain walls inter- 
polating between two adjacent vacua. Above a first critical mass m*, these BPS walls 
combine and form one non-BPS wall. If the mass is higher than a second critical value 
m**, there is no domain wall left; the two adjacent vacua can only be connected by a two 
wall configuration passing through the chiral symmetry conserving vacuum. In Ref.[|3]] 
it was found that for any theory with Nf/N c < 1/2 there is one BPS wall connecting 
adjacent chiral symmetry breaking vacua both in the small m limit and in the large m 
limit. 

The purpose of this work is to provide a unified description of the domain wall spec- 
trum of the VYT effective Lagrangian for any value of N c and Nf. One key observation 
we make is that, with a canonical Kahler potential, the static solutions to the equations of 
motion take the form 

(f) = <f>(N c z,m,Nf/N c ), 

X = x(N c z } m } N f /N c ). (1) 

Here <p and \ are me scalar components of the superfields $ and X respectively. Theories 
with different values of N c but identical ratios of Nf /N c have domain wall solutions that 
differ only by a scale factor N c . Therefore, only theories with mutually prime values of 
N c and Nf need to be considered. Theories for which N c and Nf have common divisors 
are related to theories with lower values of N c and Nf. 

We also make use of the fact that BPS states lie in shortened multiplets of the cen- 
trally extended supersymmetry algebra. Only two shortened multiplets can combine to 
form one regular multiplet, and therefore the number of BPS saturated domain walls 
can only change by a multiple of two when the mass parameter is continuously var- 
ied [^p. For small values of m, we show that if Nf/N c < 1/2 there is one BPS wall, 
and if 1/2 < Nf/N c < 1 there are two BPS walls. It then follows that the number of 
BPS walls is odd for any value of m if Nf/N c < 1/2, and even for any value of m 
if 1/2 < Nf/N c < 1. Note that this is fully consistent with the previous results in 
Refs.[g,|l. 



1 



Even though the single BPS wall at small m for Nf/N c < 1/2 can not disappear 
when the mass is increased, one can not exclude the possibility of the appearance and 
subsequent disappearance of a new pair of BPS walls at intermediate values of m. We 
therefore perform a numerical study of this class of theories, and find that such an exotic 
event does not occur. 

For theories with 1/2 < Nj/N c < 1, the fact that the number of BPS wall is even 
does not block the disappearance of the pair of degenerate BPS walls that exist at small 
m. Indeed, in the specific cases previously studied [Q] this happens at the mass m*. We 
also perform a numerical study of the theories in this class, and find that the prototype 
behavior found in Ref.[[Zp is valid in general. These theories all have a critical mass m* 
at which the two BPS walls combine to form a non-BPS wall, and a second critical mass 
m B , where the non-BPS wall disappears. It is clear from the general form of the static 
solutions in Eq.([T]) that the values of m* and m** only depend on Nf and N c through the 
ratio N f /N c . 

The remainder of this paper is organized as follows. In Section |T] we introduce the 
VYT Lagrangian as an effective description of SQCD. Next, in Section [| we discuss 
domain walls in the VYT model, followed by the presentation of our numerical results in 
Section |3| and |]. We first show that in the theories with Nj/N c < 1/2 there is one BPS 
wall for all values of m. Then, for the theories with 1/2 < Nf/N c < 1, we determine 
and 771** as a function of Nf/N c . In the final section we discuss how our results depend 
on some of the assumptions we made, and to what extend the results we derived in the 
context of the VYT effective Lagrangian apply to SQCD. 



1 Super symmetric QCD and the VYT effective Lagran- 
gian 

The fundamental Lagrangian describing SQCD, extensively discussed in for example the 
review [§], takes the form 

C = (J-^ tr J d 2 9 W 2 + h.c.J + d4d Q t] z V Qi + \Jd A e e~ v Q^j 



d 2 9 QiQj + h.c. 1 . (2) 



Here color indices have been suppressed, and i = 1, . . . , Nf is a flavor index. The vector 
superfield V contains the gauge bosons and gauginos, and Qi and Qi are the quark chiral 
superfields transforming respectively under the N c and N c representations of the SU(N C ) 
gauge group. In addition, g is the gauge coupling, and m u is the quark mass matrix. If 
Nf < N c — 1 then the theory has a dynamically generated superpotential. In the absence 
of a quark mass, the model has a SU(N f ) <g> SU(N f ) ® U(l) v <g> U(1) A <g> U(l) R global 
symmetry at the classical level. The U(\)r symmetry is broken down to the discrete 
symmetry Z, 2 n by an anomaly. It is further broken down spontaneously to Z 2 when a 
gaugino condensate (tr AA) is formed. 
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The VYT effective Lagrangian proposed in Ref.[jT]] and modified in Ref.[^] to properly 
realize the discrete anomaly free Z 2 n symmetry, provides a partial description of SQCD 
with Nf < N c at low energy in terms of the composite chiral superfields 



32vr 2 



tr W z 



Mij = 2Q t Q 



(3) 



representing the glueball and matter chiral superfields, respectively. 

The superpotential of the model, which is constructed to reproduce the anomalous 
Ward identities of SQCD, takes the form 



W = -$ 3 
3 



log 



( A ss.-s, ~(^-N f ) 



- tr (mM) . 



(4) 



The VYT effective Lagrangian is not an effective theory in the usual Wilsonian sense. It 
does not describe all relevant low energy degrees of freedom, but it reproduces Green's 
functions involving only the $ and the M fields [§]. 

We take the quark mass matrix to be diagonal and degenerate, m ,J = mb l K We then 
assume that for the purpose of studying the vacuum structure and domain walls, we can 
take Mij = X 2 5ij. In terms of X, the superpotential reads 



W 



log 



/$3(N c -N f ) X 2N f 
\ A3N C -N f 



- {K - N f ) 



m 



N f X\ 



(5) 



The Kahler potential is not uniquely determined by symmetry considerations. For now, 
we take it to have the canonical form 



K. = + XX. 



(6) 



In Section |5| we will comment on how our results may change with a different form of the 
Kahler potential. The kinetic terms for the scalar components are 



C kin = d^d ll <j ) * + d^xd,x\ 
and the scalar potential is (A = 1 from here on) 

V(0, X )=4|0 2 log fyWo-N,)^,] | 2 + N 2 



mx 



3X 



(7) 



(8) 



This potential has iV c chiral symmetry breaking vacua, corresponding to the physical 
expectation values 



N f /N c 



^2-Kil/Nc 



(x 2 h = l^- 



4 ^(N -N f )/N c 



3m / 



(9) 
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where I = 0, .., iV c — 1. In this paper we consider domain walls that interpolate between 
such vacua. In order to ensure that the two connected vacua lie on the same sheet of the 
logarithm in the scalar potential (Q), we write the vacuum expectation values of <fi and x 
as 



N f /3N C 



3m \ 



, A K(N C -N f )/2N C 

(X)k = {x)^ mk{N - Nf)IN % ^h=(^) . (10) 

with k = 0, .., N c — 1. This ensures that the phase of the argument of the logarithm is 
equal to zero in the vacuum. Note that not all vacua can be written in the form of Eq.([T0|). 
For example, the vacuum with I = 1 in Eq.@ in the case N c = 4 and Nf = 2 does not lie 
on the same sheet as the vacuum with 1 = 0. The latter vacuum is obtained both for k = 
or 2 in Eq.([T0|). These two vacua can therefore not be connected by a domain wall, unless 
it crosses the branch cut. The SQCD gaugino condensates in the vacuum labeled by k is 

(tr AA) fc = ^(0 3 ) fe = (^y^ e -^ f/Nc (U) 

Apart from the chiral symmetry breaking vacua, the VYT model also has a symmetric 
vacuum with = and \ — 0> corresponding to vanishing gaugino condensate The 
existence of a corresponding vacuum in SQCD is very controversial. 

In terms of the physical field S = $ 3 , it is clear that the chirally symmetric vacuum in 
the VYT Lagrangian arises from a singularity in the metric that is derived from the Kahler 
potential (||), whereas the chiral symmetry breaking vacua are associated with extrema of 
the superpotential (|5]). As noted before, the Kahler potential is not completely determined 
by symmetry considerations, in contrast to the superpotential. Results that depend on the 
particular details of the Kahler potential are therefore suspect. 

In favor of the chirally symmetric vacuum, it was argued in Ref.[§] that the existence 
of such a vacuum in combination with a vacuum averaging hypothesis could resolve a 
long standing inconsistency concerning gluino condensates calculated using the strong 
coupling instanton approach on the one side, and the weak coupling instanton approach 
on the other. However, in Ref.[|7p it was shown that a similar inconsistency persists in su- 
persymmetric QCD with light matter in the adjoint representation, even with the vacuum 
averaging hypothesis. As this theory is equivalent to softly broken J\f = 2 supersymmet- 
ric Yang-Mills theory (SYM), there is no uncertainty about its vacua if the mass of the 
adjoint is small. The result therefore casts doubt on the motivation for the chiral symmetry 
breaking vacuum in supersymmetric QCD. 

In Ref.[|]] it was stated that in the chiral symmetry conserving vacuum of the VYT 
Lagrangian certain discrete anomaly matching condition can not be satisfied. In Ref.[^] 
on the contrary it was argued that such considerations do not constrain the existence of a 
symmetric vacuum. 
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Another piece of evidence pointing against the existence of a chirally symmetric vac- 
uum in SQCD stems from deformed J\f = 2 supersymmetric theories. J\f = 2 SYM 
theory with a small mass for the adjoint has only two vacua. TV = 1 SYM is obtained 
from this theory in the limit that the soft breaking mass goes to infinity. In Ref.[[K| [11]] 
it was thus argued that J\f = 1 SYM also has only two vacua, barring a phase transition. 
This argument indicates that there is no room for a third, chirally symmetric vacuum. A 
similar argument concerning SQCD with one flavor can be made by deforming J\f = 2 
supersymmetric QCD with one hyper multiplet []T2]]. 

We conclude that it is dangerous to extend results obtained from the VYT effective 
Lagrangian to SQCD if these results depend on the existence of the chirally symmet- 
ric vacuum. In this article we only consider domain walls interpolating between chiral 
symmetry breaking vacua. 



2 Domain walls in the VYT effective Lagrangian 

Domain walls are static field configurations that interpolate between two vacua. They are 
static solutions to the Euler-Lagrange equations 

_ _dV_ cj\ _ _9V_ 

dt 2 dz 2 d<p* ' dt 2 dz 2 dx* 

It can be checked that static solutions to these equations take the form of Eq.(jl]). We will 
refer to the theory with N c colors and Nf flavors as (N c , Nf). It is clear that theories with 
the same ratio of Nf/N c have domain walls that only differ by a scale factor. 

BPS saturated domain walls [ p"5| ] in the VYT model satisfy the first order equations 

dz dip* ' dz dx* 
Solutions to the BPS equation have an "integral of motion" 

I = lm{e~ iS W}, (14) 

which does not depend on the space coordinate z. 

We consider domain walls interpolating between adjacent vacua. For definiteness we 
take the vacua for k = and k = 1. Note that these vacua are only neighboring in terms 
of the phase of the gaugino condensate if Nf = 1 or Nf = N c — 1. Evaluation of the 
constant of motion / in these two vacua provides the value of the phase 5 in the BPS 
equation, 

where the choice of the sign determines which of the two vacua is located at z = +oo, 
and which at z = — oo. The BPS bound on the energy of the interpolating domain wall is 
given by 

/4\ l ~ Nf/Nc , f Nt\ 
E BPS = 2N C ( | J sin ^jfj . (16) 
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Figure 1: The solution to Eq.© for the (5,1) (dotted), (5,2) (dashed), (5,3) (dotted- 
dashed) and (5,4) (continuous) theory. Open triangles indicate the m = domain wall so- 
lution^), while open circles represent the domain wall solution corresponding to m = m*. 



We parametrize the domain wall profiles in terms of four real functions by means of 



(0) o i2(z)e*M 



(17) 
(18) 



If the domain wall profiles have reflection symmetry with respect to their center, then 
the integral of motion ( |Tjj ) evaluated at the center (z = 0) yields the following relation 
between the values of R = R(0) and p = p(0) 



^/ 2 



COS ( TT^- j . (19) 



Note that iV c and Nf enter in Eq.(|T9]) only through the ratio Nf/N c , and there is no 
dependence on m. A BPS saturated domain wall in the theory (N c , Nf) correspond to a 
point on a curve defined by Eq.(|T9"|) in the pi versus R^ plane. There is one such curve 
for every value of Nf/N c . In order to determine exactly where on the curve the point is 
located, the BPS equation needs to be solved explicitly. When the mass m is varied and 
the domain wall remains BPS saturated, the point representing the wall will move along 
its curve. 

We now consider the small mass limit m <C 1, as previously discussed in Refs.[|T?[ Bp. 
In this limit the field $ can be integrated out by imposing the condition d\V/d(fi = 0. In 
the resulting theory, which only contains the field X, there is one BPS wall interpolating 
between the k = and k = 1 vacua. This wall corresponds to the maximum[] of the curve 
Pq versus Rq defined by Eq.(|T9"|). However, for certain values of iV c and Nf there is a 



'it is amusing to observe that the maxima of the curves of the (N c , Nf) and (N c , N c — Nf) theories are 
related by reflection in the line p$ = Rq. 



6 



second possible BPS wall in the small mass limit for which = 0. This wall corresponds 
to the point 



As po is a modulus and therefore larger than or equal to zero, the second BPS wall is only 
present when N f /N c > 1/2. 

We thus see that for small m there is one BPS saturated wall if Nf/N c < 1/2, and 
there are two BPS walls if 1/2 < N f /N c < 1. Therefore, if N f /N c < 1/2 the number of 
BPS walls must be odd for any value of the mass parameter m. In particular, the wall at 
small m can not just disappear. Similarly, if 1/2 < Nf/N c < 1, the number of BPS walls 
must be even. This allows the possibility of the simultaneous disappearance of the pair of 
BPS domain walls that exists at small m. 

We illustrate this analysis in Fig. [IJfor the theories with N c = 5 and Nf = 1, 2, 3 and 
4. We draw the four curves given by Eq.(|l9t) in the p\ versus Rq plane, and indicate the 
location of the BPS walls at small m on the curves with a triangle. For Nf = 1 and 2, 
there is only one triangle on the curve. When the mass m is increased, these triangles 
start to move along the curve. How they move is a dynamical question which can only 
be answered by explicitly solving the BPS equations. It is clear though that they have 
to remain on their curve, they have nowhere else to go. In principle, there is no a priori 
argument that precludes the appearance of a new pair of BPS walls. As we will show in 
Section this does not occur. For Nf = 3 and 4 there are two triangles on the curves. As 
we will show in Section |], when m is increased, the two triangles start to move towards 
each other, and at the critical value m* they meet at the point on the curve indicated by a 
circle. They then combine to form a non-BPS wall. 



The existence of one BPS wall at small m for theories with Nf/N c < 1/2, in combination 
with the fact that the number of BPS walls modulo two is conserved, implies that the 
number of BPS walls in this class of theories is odd for any value of m. In order to 
determine whether any new pairs of BPS walls or additional non-BPS walls appear when 
m is increased from zero, we performed a numerical simulation of all the theories in this 
class with N c < 8. 

We simulated the second order equations of motion on a lattice with a forward predict- 
ing algorithm. To that purpose, we transformed the Euler-Lagrange differential Eqs.([T2|) 
to finite difference equations, with discretized space and time coordinates. The scalar 
field <p and x were kept fixed at their vacuum expectation values at the two ends of the 
lattice. The size of the lattice is chosen to be much larger than the width of the walls 
and, at the same time, the lattice spacing is chosen much smaller than the width. An 
initial configuration for and x was constructed so that it interpolates between the two 
vacuum expectation values. This initial configuration does not satisfy the static equations 
of motion and therefore evolves in time. We are interested in a minimum energy con- 



R 





(20) 



3 The case N f /N c < 1/2 



7 



0.5 



-0.1 0.1 0.2 




2 . 5 

2 

1.5 

1 

0.5 





0.1 0.2 

Z 



Figure 2: BPS domain walls for the (5,1 ) theory. Masses are m — 10 (dotted), 20 (dashed) 
and 50 (solid) 



figuration. The system can only reach such a minimum energy configuration from the 
initial configuration if it can dissipate energy. We therefore added a friction term to the 
equations of motion. The coefficient of friction is adjusted so that the system relaxes to 
the minimum energy configuration in as short a time as possible. If the coefficient is too 
large, the system does not converge, and if it is too small, the systems keeps oscillating. 
Of course, as the friction terms are proportional to the time derivatives of the fields, the 
static configurations of the equations of motion with and without the friction term are the 
same. Once the system had come to rest, we calculated the energy of the final configu- 
ration and compared this to the BPS bound. We also explicitly checked that domain wall 
solutions we thus obtained did not cross the branch cut of the logarithm appearing in the 
scalar potential (||). This check was performed by plotting the argument of the logarithm, 
^3(N c -N f )^2N f ^ m com pi ex plane. As the fields on the domain wall vary from one 
vacuum to the other, this quantity traces a loop in the complex plane. We then verified 
that this loop does not cross the negative real axis (where the branch cut is located within 
our conventions). 

We repeated this procedure for each theory, starting at low values of m and following 
the walls as m is increased. In order to see if more than one wall exists, we also repeated 
the procedure with various initial conditions. We found that for the theories in the class 
Nf/N c < 1/2 there is only one BPS saturated domain wall for any value of m. In Fig. ^ 
we show some representative examples of BPS domain walls for the (5,1 ) theory and 
various values of m. Following the domain wall in this theory in Fig. |T|, it starts on top of 
its curve at m = 0, and then moves along the curve towards smaller values of when m 
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Figure 3: The difference between the energy E of the domain wall configuration and the 
BPS bound -Ebps a s a function of the mass m (dots) for the (5,3) theory. The dashed line 
indicates the difference _E 2W — E BPS . 



is increased. 



If 1/2 < Nf/N c < 1 there exist two BPS domain walls in the small m limit. The 
argument that the number of BPS walls can only change by a multiple of two allows for 
the possible existence of two critical masses, m* and m**, as found in Ref.[@] for the 
theories (N c , N c —1) and N c = 2, 3 and 4. At the mass m*, the two BPS walls combine to 
form one non-BPS wall, and at the mass m** the non-BPS wall disappears. We perform 
a numerical simulation of all theories in the class 1/2 < Nf/N c < 1 with N c < 8, and 
find the same sequence in each theory. The theories we simulated provide a sample of 
values of Nf/N c that covers the whole allowed range of 1/2 < Nf/N c < 1. We therefore 
conclude that the theories in this class, not just the ones with N c < 8, behave qualitatively 
in the same way. Of course, the actual values of m* and m** depend on Nf/N c . 

In Fig. [3] we show how the difference between the energy of the domain wall and the 
BPS bound depends on the mass m in a representative case, the (5,3) theory. We have 
indicated the masses and m*#, and we also plot the sum of the energy of the two BPS 
walls that connect the two vacua for k = and k — 1 through the intermediate vacuum at 
= and x = 0, 



Below m*, the energy of each of the two possible walls is equal to the BPS bound. When 
the mass is increased above m„ the energy of the non-BPS wall starts to increase and 
tends towards E 2 w- In feet, at a mass below m** the energy of the wall exceeds E 2 w- 
However, at that point there is still a sphaleron barrier, and the wall is meta-stable. The 



4 The case 1/2 < N f /N c < 1 
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Figure 4: Logarithm of the critical mass m* versus the ratio Nf/N c 



height of the barrier decreases as the mass approaches m**, and when the mass is equal 
to m**, the non-BPS wall ceases to exist. 

Let us pause briefly to account for the kink in between and m»* appearing in Fig. |3|. 
At this value of the mass parameter m the wall changes rapidly from a shape resembling 
the BPS configuration, to a shape more similar to a two wall configuration approaching 
the vacuum in the middle. 

In order to find m* numerically, we start with a small value of the mass parameter m 
for which the energy of the domain wall configuration does not differ significantly from 
the BPS bound. The initial configuration determines which of the two possible walls we 
obtain; this does not matter. Then we gradually increase m and determine m* as the value 
of m for which the energy of the wall differs significantly from the BPS bound. 

In Table [I] we present the numerically determined values of the critical mass m* for 
the (N c , Nf) theories with iV c < 8. We also plot the logarithm of as a function of the 
ratio Nf/N c in Fig. |]. We are not able to determine numerically for Nf/N c > 4/5 
because its value becomes too small. However, from Fig. |] it appears that this critical 
mass goes to zero in the limit Nf/N c — > 1. 

In order to determine the critical mass m„, we use the following procedure. For a 
fixed value of m we calculate the minimum energy E for field configurations interpolating 
between the given vacua subject to a constrained value of R . If m is just above m*, a plot 
of E versus R shows a minimum and a maximum, corresponding to a non-BPS domain 
wall and a sphaleron, respectively. As m increases the minimum and the maximum move 
towards each other and become less pronounced. We determine m = m** as the value 
of m were the maximum and the minimum coincide. For higher values of m, the energy 
E increases monotonously with R , which means that the non-BPS wall has disappeared 
and the minimum energy configuration is formed by two walls connecting through the 
chirally symmetric vacuum. 

In Fig. |5] we illustrate this procedure for the (5,4) theory. The numerically determined 
values of the critical mass for various (N c , Nf) theories are presented in Table ^|. In 
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Figure 5: The minimum energy E of interpolating field configurations subject to a con- 
strained value of Rq in the (5,4) theory. The arrows indicate the position of the non-BPS 
domain wall (solid) and of the sphaleron (dashed). As m increases, the minimum and the 
maximum move towards each other while the height of the barrier decreases. At m„ the 
wall and the sphaleron coincide and then they disappear. 



Fig. H we plot in** versus the ratio Nf/N c . It can be seen from this figure that in the limit 
Nf/N c — * 1, m** tends to the constant value m** rs 1.83. 



5 Discussion 

We derived our results in the context of the VYT effective Lagrangian. To what extend 
these results are valid for SQCD therefore depends on how well the VYT effective La- 
grangian describes the vacuum structure of SQCD. An interesting discussion of this issue 
is given in Ref.[§]. 

In addition, we assumed that the Kahler potential for the fields $ and X has the 
canonical form @. The tension of the BPS walls does not depend on the Kahler potential, 
but the shape of the domain walls and quantities such as the critical masses m* and m** 
do depend on its details. It is instructive to see what happens in the case of a simple 
deformation of the Kahler potential, 



K = a ($$ + PXX) 



(22) 



After rescaling z' = z/a, X' = ^f]3X, m! = m//3 and A' = Aj3 Nf/( - 3Nc ~ Nf \ the static 
equation of motions of the scalar components are cast back into exactly the same form 
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Figure 6: The critical mass to** versus the ratio Nf /N c . It can be seen that m** approaches 
a constant value in the limit Nf/N c — > 1. Extrapolating from our calculated points we 
find this limiting value to be to** rs 1.83. 



they have when a = /3 = 1, and our results apply to the primed quantities. Note that 
in principle both a and (3 can depend on Nf and N c , and not necessarily only through 
the ratio Nf/N c . Even when more complicated deformations of the Kahler potential are 
introduced, it is expected that the qualitative picture remains intact, at least as long as 
the singularity structure is not changed. In particular, for any value of to, the number of 
BPS walls will remain odd for N f /N c < 1/2, and even for 1/2 < N f /N c < 1 under 
such deformations. Given the suspicion surrounding the chirally symmetric vacuum, it 
would be interesting to see what happens if our analysis were repeated with a Kahler 
potential that does not contain a singularity. The fate of the lower branch for theories with 
Nf/N c > 1/2 is uncertain in that case [ |H3[ [TTp , and our conclusions could potentially be 
changed. 

Finally, the difference in the domain wall spectrum of the theories with Nf/N c < 1/2 
and the theories with 1/2 < Nf/N c < 1 poses the following question. What happens 
when the mass of one of the quark flavors in a theory with N c = 3 and Nf = 2 is moved 
to infinity? The theory with two degenerate masses has an even number of BPS walls. If 
the mass of one of the quarks is taken to infinity while the other mass is kept finite, the 
heavy quark can be integrated out and the theory with N c = 3 and Nf = 1 is obtained. 
However, this latter theory has only one BPS wall. This change in the number of BPS 
walls seems to be inconsistent with the argument that the number of BPS walls can only 
change by a multiple of two when a mass parameter is varied continuously. Obviously, 
one of the assumptions of the argument is violated. We intend to come back to this issue 
in a future work. 
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Table 1: Numerically determined values of the critical mass m* for different (N c ,Nf) 
theories. The symbol 'x' indicates that the corresponding theory has one BPS wall for 
any value of the mass parameter m and therefore there is no m*. The question mark 
indicates that the value of m* is too small to be determined by our method. 
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Table 2: Numerically determined values of the critical mass m**. The symbol 'x' indi- 
cates that the corresponding theory has one BPS wall for any value of the mass parameter 
m and therefore there is no m**. 
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